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Abstract: A disordered version of the one dimensional asymmetric exclusion model where 
the particle hopping rates are quenched random variables is studied. The steady state 
is solved exactly by use of a matrix product. It is shown how the phenomenon of Bose 
condensation whereby a finite fraction of the empty sites are condensed in front of the 
slowest particle may occur. Above a critical density of particles a phase transition occurs 
out of the low density phase (Bose condensate) to a high density phase. An exponent 
describing the decrease of the steady state velocity as the density of particles goes above the 
critical value is calculated analytically and shown to depend on the distribution of hopping 
rates. The relation to traffic flow models is discussed. 
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Introduction: Nonequilibrium systems may be defined as those evolving to some dynami- 
cal rule which a priori does not respect detailed balance with respect to an energy function. 
The interest in such systems is that they describe a broader class of phenomena, for ex- 
ample driven diffusive systems [p]] and growth processes 0, than can be handled within 
the traditional framework of equilibrium statistical mechanics. The asymmetric exclusion 
model is a prototypical nonequilibrium model. In various guises it has been used to model 
interface growth, polymer dynamics and traffic flow [JTJ, ^|, ^, ^J] . It comprises classical parti- 
cles stochastically making nearest neighbour hops with hard core exclusion and a preferred 
direction imposed. Through iteration of the dynamical rules a nonequilibrium steady state 
is attained. Even in one dimension steady states of asymmetric exclusion processes have 
revealed such diverse phenomena as long range correlations, shock formation, boundary- 
induced phase transitions and symmetry breaking [|, §, [7], || ||. Significant progress has 
been made in understanding the one dimensional model through the exact solution of steady 
states via a matrix technique |§. Recently this technique has been extended from continuous 



time (random sequential) dynamics to a parallel updating scheme JlO 



In the present work the effects of disorder in the model are considered. Specifically, each 
particle has its own intrinsic rates for hopping forwards and backwards which are taken 
as random variables. Previously, single inhomogeneities such as a fixed blockage [|ll[] or a 
moving blockage |12j were studied and phase transitions found. In another context a repton 
model with two allowed hopping rates was solved and a mapping to a Bose gas noted [||]. 
Here we employ a similar mapping and report that in the fully disordered model studied, 
Bose-Einstein condensation occurs. In the limit of a single inhomogeneity the condensation 



transition corresponds to a limit of [12 



To understand, at a descriptive level, the Bose condensation found here it is useful to 
visualise a caricature of traffic flow. Consider a one lane road with cars moving along, each 
at its own chosen speed, and no overtaking allowed. After some time a tailback will be 
formed behind the slowest car with all the other cars' speeds reduced to that of the slowest. 
However, if the density of the cars on the road is increased, for example, by cars merging 
with the road, a congested situation may be reached wherein the speed of all cars, regardless 
of their own preferred speeds, is limited by the high density of traffic. To see the analogy 
with Bose condensation one relates the distribution of available road between cars to the 
distribution of Bosons within the states of an ideal Bose gas. In the low density case a finite 
fraction of the available open road is condensed in front of the slowest car whereas in the 
high density situation the available road is spread thinly between all cars. 

In the following the analogy between a phase transition in the model considered and 
Bose condensation will be made exact through the exact calculation of the steady state of 
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the model. This is performed with the aid of a matrix product ansatz. A critical exponent 
describing how the velocity of a particle decreases as the density goes above the critical 
value is calculated exactly and related to the distribution of car velocities. 

Model Definition and Construction of Steady State: The model studied here is 
defined as follows. Consider M particles labelled fi = 1, . . . , M hopping on a lattice of size 
iV with periodic boundary conditions. A particle // at site % attempts hops forward to site 
2 + 1 with rate p^ and hops backward to site i — 1 with rate g M . Due to exclusion, hops only 
succeed if the target site is empty. Since particles cannot overtake the sequence of particles 
is preserved. We summarise the allowed exchanges at a pair of sites i, i + 1 as 



/U — > /J, with rate p^ ; fi — > /i with rate g M 



(1) 



In order to obtain the steady state of this model a matrix product || was used. The ansatz 
is to represent the weight (unnormalised probability) of a configuration by the trace of 
a product of matrices D^, E, one matrix for each site; we write for site % if particle 
fi is at site % in the configuration and E if a hole is at site %. For example, the weight 
/jv(wi,ri2, • • • ,tim) °f a configuration of TV sites comprising particle 1 followed by n\ holes; 
particle 2 followed by n 2 holes and so on, is written as 



f N {m, n 2 , ... , n M ) = Tr [D\ E ni D 2 E n2 . . . D M E 



(2) 



The weights are related to probabilities P/v({n M }) via a normalisation Z NjM defined by 
^iv({7fy}) = /jv({%})/^jv,m- We see from @ that in matrix form 



(N-M) 



Tr 



ni,n 2 ...n M 



' M 

n D v E n - 



(3) 



(so that the total number of holes is iV — M). A sufficient condition for (§) to indeed give 
the steady state is 

p^D^E-q^ED^ = D^ for /i = 1 . . . M (4) 

Details of the proof will be given elsewhere |I3| |. Due to the fact that the particles cannot 
overtake each other, the velocity v (the steady state average of the rate of forward hops 
minus the rate of backward hops) is the same for each particle fi. Taking /i = 1 



^ 5 E m «m, (^v-A/-i) Tr 
Zn-i,m/Zn,m 



M 



{p 1 D 1 E - qiEDi}E ni \{ D v E n » 



u=2 



/z. 



N.M 



(5) 
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where we have used (f|) and The existence of a non-trivial representation of @ 
provides a check that the algebra @ is consistent. However, rather than using a represen- 
tation to perform calculations we can proceed directly using Q). We begin by using (|j) for 
/i = M in (H) to write 



—f N (rii, . . .,n M -i + l,n M - 1) = — /v-i(^i, • • .,n M -i,n M - 1). 
Pm Pa/ 



/jv(ni, . . . ,n M ) 
The procedure is continued using, in sequence 



M-l 



n „ 

i=Q Pm- 



/iv(ni, . . . ,n M ) 



with /i = M-l,M-2...1to obtain 

9j 



M-l -, M 

^ _L_ j| « 

i=0 PAf-t j=Af+l-i ^ 



/jv-i(ni, . . . ,n A/ - 1) 



The effect of this manipulation has been to commute a hole initially in front of the particle 
M backwards one full turn around the ring. The result is that the weight of a configuration 
of size N is expressed as a multiple of the weight of a configuration of size N — 1 with one 
hole less in front of particle M. Repeating the commutation procedure for a hole initially 
in front of particles M — 1, M — 2, ... 1 we deduce that, to within a multiplicative constant, 
the weights are 

-l 



M 



/jv(wi, 



n M 



Way where 



Ai=l 



M-l 



E 



n 



9j 



o j= M +i-i Pi 



M 



i-n- 

k =i Pk 



(6) 



An immediate consequence of the form (|6]) is a mapping to an ideal Bose gas. The 
N — M holes of the exclusion process are viewed as Bosons which may reside in M states 
with energies determined by the hopping rates of the M particles i.e. exp(—(3E^) = g^. 
A hole in the fi th Bose state is in front of particle \x (and behind particle fi + 1). The 
normalisation Zjv,m is equivalent to the canonical partition function of a Bose gas. The 
velocity is a ratio of partition functions of different system sizes (^) and corresponds to a 
fugacity. 

One sees from that the energy levels are, in general, complicated functions of the 
hopping rates. However, a simplification occurs if the particles are restricted to hop forward 
only i.e. = V/i. In that case the algebra (Q) reduces to p^D^E = and it is easy to 
see that = For simplicity this case will be discussed in the following. 

An illustrative example of Bose condensation: To illustrate the phenomenon of Bose 
condensation in this model consider first the simplest case of an inhomogeneous system: 
particle 1 has hopping rates Pi = p , q\ = while the other M — 1 particles have hopping 
rates p^ = 1 — for fi > 1. This situation is, in fact, a limit of the problem of a 
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single moving defect studied recently fTJJ. The normalisation Zn,m is easy to calculate 



combinatorically since gi = p 1 and for \i > 1, — 1 thus 

ZN,M^{ N M n ^~ 2 ) P - n \ (7) 
rai=0 \ / 

yielding an exact expression for the velocity through (||). In the thermodynamic limit 

N -> oo with M = pN , (8) 

where the density p is held fixed, the sum (0) is dominated by n\ ~ O(N) for p < 1 — p and 
ni ~ 0(1) for p > 1 — p and it can be shown that 

for p < 1 -p Z NM ~ p-( N ~W (i _ p)-(M-i) and y _^ p (g) 
for p > 1 - p Zjv,m - f I P 2 vl (P - 1 + P) and u — >• 1 - p (10) 



Two phases are apparent: a low density phase p < 1 — p where the velocity is equal to p, 
hence determined by the slowest particle; a high density phase p > 1 — p where the velocity 
is 1 — p, hence determined by the density. The low density phase correpsonds to a Bose 
condensate since in the steady state (ni)/N (the fraction of the holes in front of the slowest 
particle) is non-zero. Equivalently, the density of holes at the site immediately in front 
of the slowest particle is 1. In the high density phase the density of holes at this site is 
(1 — p)jp. At the transition (p c = 1 — p) one may think of the tailback behind the slow 
particle reaching all the way round the ring. Above p c = 1 — p the velocity decreases from 
p according to v = p — (p — p c ) a where a — 1. 

The Disordered Case: We now turn to the general problem where each particle's hopping 
rate p^ is a random variable. To keep the analysis simple we assume particle 1 has the lowest 
hopping rate which we take to be p\ — 1 while the other particles have hopping rates drawn 
from a distribution V(p). We require that V(p) vanish at p = 1 (so that particle 1 is indeed 
the slowest) and as p — > oo, V(p) vanish fast enough that the distribution be normalisable. 
Here a simple choice related to the gamma distribution is considered 

V(p) = (p- l) 7 exp(-(p - 1)) [r( 7 + IF 1 (11) 



where IX7) is the usual factorial function. Other distributions may be studied [[TJ 



To analyse the disordered case one wishes to obtain properties given by a typical real- 
isation of the disorder (here the particle hopping rates), for example the typical velocity. 
Usually in the theory of disordered systems one expects quantities such as v to self-average 
but quantities exponentially large in the system size, such as Z NM , not to. That is, one 
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expects as iV — > oo, v — >• t> 7^ Zn-i,m/Zn,m where the bar indicates an average over the 
particle hopping rates. To circumvent this difficulty it is convenient to work with the grand 



canonical ensemble [15] where in the thermodynamic limit, we have 



where N -M = V — — r (12) 



In the thermodynamic limit ([§]) the fraction of particles with hopping rates between p and 
p + dp converges to V(p)dp. Therefore (|12"D may be replaced by 

(ti . \ poo V 

1 - p = pl( z ) + 1-±L where Hz) = dpVip) (13) 

N Ji p — z 

and (rii) = z/(l — z) denotes the steady state average of n\ (the number of holes in the 
Bose ground state). Clearly, V(p) is analogous to the density of states of the Bose gas. The 
condition for Bose condensation [|l^] translates to 

1(1) < (1 - p)/p . (14) 

Recall from the theory of Bose condensation |15] that since we must have z < 1, ([14]) implies 



(ni)/N is non-zero and z = 1. If (|14]) does not hold, then in Cp~3|) , (ni)/N = and one 
solves for z < 1. The integral I(z), with P(p) given by fllT]), may be expanded in powers of 
1 — z via a Mellin transform |T3tl to obtain 



for 0<7<1 I(z) = I/7 + IX-7) (1 - Z y + ... (15) 
for 7 > 1 I(z) = I/7 - 1/(7 - 1) (1 - z) + . . . (16) 

Hence, by ( |I4"D one has p c = 7/(1 + 7) and just above p c , v decreases as 

for < 7 < 1 v ~ 1 - [-r(- 7 )p;?]~ 1/7 (p _ p c )V7 ( 17 ) 
for 7 > 1 v ~ 1 - (7 - l)p^ 2 (p - p c ) (18) 

For p < p c we have Bose condensation: the fraction of holes in front of the slowest particle 
is (ni) / (JV — M) = 1 — 7(l)p/ (1 — p). At the transition this fraction vanishes and above the 
transition the velocity decreases from 1 as 1 — v ~ (p — p c ) a , where = 1/7 for 7 < 1 and 
a = 1 for 7 > 1. For 7 > 1 the same a as for a single slow particle is found. 

Discussion: In this work it has been shown that the steady state velocity in a disordered 
exclusion model undergoes a phase transition exactly equivalent to Bose condensation. In 
the low density phase (Bose condensate) the velocity is equal to the hopping rate of the 
slowest particle. Above a critical density a high density phase is entered and the velocity 
decreases according to an exponent which has been calculated exactly here for the totally 
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asymmetric case and a specific distribution of hopping rates flTTp . The exponent depends on 
the distribution of slow hopping rates which may be considered analogous to the density of 
states of the Bose gas. Other distributions of hopping rates yield a similar behaviour ||13|| . 

In the latter part of this letter the scope was restricted to a totally asymmetric exclusion 
process where only forward hops are allowed. This implied that in the Bose gas the ener- 
gies were independently distributed, reminiscent of a Random Energy Model of disordered 
systems (TIJ . It would be of interest to investigate more fully the partially asymmetric case 
where backwards hops (g M 7^ 0) are allowed. As is seen from (||), in that case the energies of 
states in the Bose gas are correlated. In the present work the ground state, i.e. the slowest 
particle, was taken to be non-degenerate. In the case of high degeneracy of the slowest 
particle, such as only two hopping rates of particle, or of a distribution V(p) peaked at 
p = 1 the phase transition is not present Jl3| . 



In the introduction the model was discussed as a caricature of traffic flow. More realistic 
but related traffic models based on asymmetric exclusion with parallel updating have been 



studied H. Interestingly, bunching of cars has been noted in these models [14 and it is 



important to determine whether this is related to the Bose condensation analysed here. 
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writing of this work and B. Derrida for helpful discussions. 
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